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^ ' ABSTRACT 

O , The atomic recombination process leads to a softening of the matter equation of state 

' as reflected by a reduced generalized adiabatic index, with accompanying heat re- 

, lease. We study the effects of this recombination softening and reheating of the cosmic 

f^) ' plasma on the ionization history, visibility function, Cold Dark Matter (CDM) transfer 

function, and the Cosmic Microwave Background (CMB) spectra. The resulting mod- 
ification of the CMB spectrum is 1/10 of WMAP's current error and is comparable to 
PLANCLCs error. Therefore, this effect should be considered when data with higher 
' accuracy are analysed. 

, Key words: Cosmology - cosmic microwave background, theory 
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1 INTRODUCTION 

l' Recent advances in the measurement of Cosmic Microwave Background (CMB) A nisotropics have attracted considerable 
interest in the field (for latest CMB observations, please refer to lBennett et alJ I^OO.'j) . website of WMAP Sate llite 1 , website of 
PLAN CK Satellite 2 , Max Tegmark's CMB experiment website 3 ; for recent reviews of CMB theory, please read lHu fc Dodelsonl 
, J2002I) and lDurrerl J200ll) ). The cosmic microwave photons are believed to be remnants of radiations produced in the hot early 
universe and carry information about the state of the universe during the recombination epoch, when the photon cross section 
, is drastically reduced due to the formation of neutral atoms. It is therefore possible to put tight constraints on cosmological 
theories using high precision CMB Anisotropics data, provided that the theoretical computation of the CMB Anisotropics 
also achieves at least the same level of precision. We need to understand the recombination process well if we want to extract 
■ information from CMB Anisotropics correctly. 

In the standard calculation of the recombination process dSeaeer et al.lll999l. 2000), the baryons are treated as ideal gas 
with an adiabatic index of 7 = 5/3. However, as the protons and electrons combine t o form neutral atoms, t he equation of state 
(EOS) of matter is softened, as reflected by a reduced generalized adiabatic index l|Chandrasekharlll939l:lMihalas fc Mihalasl 
1999), which depends on the ionization fraction. This change in the generalized adiabatic index is most drastic during 
the recombination, and the heat produced delays the completion of the process. In this paper, we study the effect of this 
recombination softening and reheating of the cosmic plasma on the ionization history, visibility function, the CMB spectra, 
and the cold dark matter (CDM) power spectra. 

We first present a brief review of the standard recombination calculation in Section [5] which is then followed by a 
discussion of the modifications of the standard equations to take into account the recombination softening and reheating in 
Section We then present our numerical results in Section 0] and summarize our discussion in Section 03 
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2 STANDARD RECOMBINATION CALCULATION 
2.1 Expansion of the universe 

The recombination process did not happen instantaneously because of finite interaction and transition time. Thus we need 
to follow the time evolution of the ionization fraction of hydrogen and helium and also the baryon temperature. The rate 
equations can be written in terms of time t, which is related to the redshift z by 

^ = -(l+z)H(z), (1) 

where H(z) is the Hubble parameter. The expansion of the universe as a whole is described by the cosmic scale factor 
a(t) = 1/(1 + z). 

The evolution of the Hubble parameter is determined by the Friedmann equation 

H{zf = H 2 [n R (i + z) 4 + n M (i + z? + Qk(i + zf + n A (i + zf (1+w) ] , (2) 

where Ho is the present value of the Hubble parameter and can be written as 100/tkmsec -1 Mpc -1 . The fi's are the density 
parameters, with subscripts R referring to the radiation contribution, M the matter contribution (including both baryon B 
and cold dark matter C), K the curvature contribution, and A the dark energy contribution, w is the EOS index of dark 
energy, which is assumed to be constant. The radiation density parameter is calculated by 

8nG{l + f v )U 

llR = o/rr x 2 ' \ 6 > 

3(Hocy 

where /„ = 3 x (7/8) x (4/ll) 4 ^ 3 ~ 0.681 is the contribution of neutrinos to the energy density of photons assuming three 
massless neutrino types, G is the gravitational constant, and U is the photon energy density given by Stefan's law U = «h,Tr. 



2.2 Ionization fractions 

We now consider the equations for the ionization fractions. Due to the higher ionization energy of helium comparing to 
hydrogen, helium recombined first. A helium atom has two electrons, and so there are two processes of recombination, namely 
Hen recombination (from He 2+ to He + ) and He I recombination ( from He + to He atom). T his Hell recombination process is 
so fast that it can be approximated by using the Saha equation JSeaeer et alJll999t 2000), at least in high-fie models. For 
1ow-S7b models, the Hell recombination is slightly slower than the Saha recombination. But since the Hell recombination 
occured too early and had negligible effect on the CMB power spectrum, it should be adequate to calculate He n recombination 
by using the Saha equation. 

The Saha equation for He n recombination is 4 

(ge - 1 - fnc)xe = (27rm e fcr M ) 3/2 - XHeU /kT M 
1 + 2/hc - x c h?nn ' { ' 

where x e = n e /nn is the ionization fraction with n e being the electron number density, nn the total number density of 
hydrogen nuclei, /h c = jihc/hh = Y p /A(l — Y p ) the ratio of the numbers of helium nuclei to hydrogen nuclei, riHe the total 
number density of helium nuclei, and Y p is the primordial He abundance by mass. Here, m e is the mass of electron, k is 
Boltzmann's constant, Tm is the matter temperature 5 , h is Planck's constant, and XHoii is the ionization potential of He 2+ . 

For H and He I recombinations, the rate equations should be used. A common me thod for the calculation of the recom- 
bination process is sometimes called the Peebles recombination JPeebleslll96l Il993h . This method considers a three- level 
atom with the ground state, first excited state and continuum. A recombination coefficient a is used in the calculation. If one 
excludes recombinations t o the ground state and assumes th e Lyman lines to be optically thick, the recombination is known 
as Case B recombination JHu et alJll995l : ISeaeer et alJl2000l) . 

Here we consider hydrogen first. The treatment for helium is essentially the same. The ionization rate is given by 

^ = C T [/fe(T M )(l - x p )e- h ^ 2 ° /kTM - nnan{T M )xeX P ] (5) 

JPeebleslll96lll99dlMa"fe Bertschingell995llSeager et alJl99sll2OO0h . Here C r is a reduction factor (to be discussed below) , 
PH2s is the Lya frequency, qh(?m) is the Case B recombination rate for hydrogen, and the total photoionization rate Ph{Tm) 
is given by 

4 This form of t he Sa ha e quation is different f rom the usual one (e.g. in lMa fc Bertschingeil ll995t) V Here we have followed the notation 
of lSeager et alJ il999ft and lSeager et all l2000ll . The left hand side is different because the ionization fractions are all defined relative to 
hydrogen. ^^^^^^^^^^ 

Matter is coupled by S-wave scattering iHannestadll200ll) . By comparing the scattering rate and the Hubble expansion rate, it can be 
proved that the assumption of a single matter temperature is valid at least down to a redshift of several hundreds. 
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(6) 



Due to the large Ly man alpha and Lyman continuum opac ities, recombination directly to the ground state will lead to 
immediate re-ionizaton |Peeblesj|l96 q : iMa fe Bertschingerlll995ri . We therefore do not include recombination directly to the 
ground state. Recombination occurs either by 2s level to ground state decay (with rate Ah), or by cosmological redshifting 
of Lyman alpha photons (with rate A Q ). Since an atom in the 2s level may be ionized before decaying to ground state, the 
recombination rate is reduced. The reduction factor C r is the ratio of the net decay rate to the sum of the decay and ionizaton 
rates from the 2s state, given by 

A tt + A H 



a 



A^+Ah+M^m) ' 



where 



A a = 



8tt da 



8irhc 



m, 



At 



3 x 13.6 eV 



(7) 



(8) 



It is a very good approximation to replace m s by (1 — xh)tih for Tm <C 10 5 K. 

Now all we need ar e the recombinatio n coefficients for H and He II. The recombination coefficient of hydrogen an is given 
by the fitting function JSeaeer et alJl200Cl) 



an = 10" 



at" 



(9) 



1 + ct d 

with a = 4.309, b = -0.6166, c = 0.6703, d = 0.5300, and t = T M /10 4 K. The recombination coefficient of He I a Ho i is fitted 



a Hoi 




i-p 




i+p- 



3 - 

m s 



(10) 



with q = iq-16 74 4 p = Q 7n ^ = 1 Q 5.ii4 K and Ta fixed arbitrarily at 3 K 

Note that in lSeager et alJ l|l999h . a fudge factor 1.14 is multiplied to Eq. E]to approximate the result of a more detailed 
calculation. The detailed calculation includes totally 609 levels. There are 300 levels of H, 200 for He I, 100 for Hell, 1 each 
for He ill, electrons, protons, and each of the five molecular or ionic H species. The purpose of our work is not to reproduce 
the 609-level result. Since we just want to investigate the effect of taking into account the ionization energy of particles, and 
using the fudge factor will give unnecessary complication of the problem, we set the fudge factor to 1. Then the whole set of 
equations reduces to the Peebles recombination. 



dx p 
~dz 

and 



We can now write down the two ionizat ion fraction rate eq uations for the proton fraction x p and the singly ionized helium 

" eager et~ai]ll999l) 

] [1 + KnAnnn(l - x p )} 



tion fraction rate cq 1 

fraction xudi = «-Hcii/n-H, as described by JSeager et alJll999ti 



[xeXpnuaa ~ /?h(1 



- h "H 2s/ kT M 



H(z)(l +z)[l + Kh (Ah + /fe)n H (l - x p )] 



dx'HcII 


[xHellXeWHaHcI — /fed (/He " 


-XHdOe-^Ho^i^/^M] 


1 + A' Hc iAHcn H (/Hc - x Hc ii)e-' 1 ' / H (! i2i P 2i s /^M 


dz 


H(z)(l + z) 


1 + Khci(A H c +/3 Hc i)n H (/H C - XHeii)e~ hVa ° 12l » 2l » /k ™ 





(11) 



(12) 



The electron fraction is x e = x v + a;Hoii- The H Lya wavelength is Ah2 P = 121.5682 nm. The H 2s — Is frequency 
is ish2s, nearly equal to c/\u2 P - The Hel 2 1 p — l x s wavelength is A Hc i2i P = 58.4334 nm. The Hel 2 1 s — l x s frequency 
is ^Hci2i s = c/60.1404 nm, and ^ H oi2i P 2i a is g ive n by ^ Hc i2i P 2i s = ^Hoi2i P ~ »^Hei2i»- The H 2s - Is two-photon rate is 
Ah = 8.22458 s _1 , and the Hel 2 1 s — l x s two-photon rate is Ah c = 51.3 s _1 . Ku = A H 2 P /[8vr//(z)] is the cosmological 
redshifting of H Lya photons, which is the reciprocal of A Q in Eq. [7| The cosmological redshifting of He I 2 x p — l 1 s photons is 
J^Hci = A^ cI 2l /[8ttH(z)]. The radiation temperature is given by Tr = To(l + z), with To being the present CMB temperature. 



2.3 Matter Temperature 

The major effects affecting the evolution of matter temperature are C ompton cooling an d adiabatic cooling. Other cooling 
and heating terms have effects at 10 _3 % level in the ionization fraction iSeaeer et alJl200oT) . and so we also ignore those terms 
in our calculation. 



There is a typo error in lSeager et alJ [l999) Eq. 4. 
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Electrons and photons are coupled by Compton scattering. When electrons and photons are nearly in thermal equilibrium, 
the rate of energy transfer is given by iPeeble sl ll968l Il993l : ISeaeer et al.ll2000h 

dE e Aa T Un e k 

—t— = (Ir — Im), (-13) 

at m e c 

where E e is the energy density of electrons and or is the Compton scattering cross section. The matter temperature therefore 
obeys: 

^ = !^^(T R -T M ), (14) 
dt 3 m e c n to t 

where n to t is the total number density of particles. 

We next consider the adiabatic cooling term. We treat the matter as an ideal gas with an adiabatic index 7 = 5/3. 
Recalling that Tm oc p 7_1 and p oc (1 + z) A (ignoring the pressure of baryonic matter in the Friedmann equation 7 ), it can be 
proved easily that the adiabatic cooling rate is described by 

^ = -3( 7 - l)H(t)T M = -2H(t)T M ■ (15) 
Therefore the evolution of the matter temperature with respect to redshift z is 

dT M _ 8a T U n e 2T M , . 

dz 3H(z)(l + z)m e c n e + n H + n Hc 1 M RJ (1 + z) ' ( ' 

By solving Eg. 1111 Eg. I12l and Eg. 1161 the ionization fractions and matter temperature can be found as a function of z. 



3 RECOMBINATION SOFTENING 



In the derivation of Eq. 1161 we have assumed that the baryonic matter is an ideal gas. It is a good approximation when the 
gas is either neutral or fully ionized. However, a partially ionized gas tends to recombine under compression, and so the EOS 
should be softer than that of an ideal gas. Here we rewrite the temperature equation and investigate the effect. 

Again we consider the recombination of H and He I only. The specific internal energy is the sum of the translation, 
excitation, and ionization energies of the particles H and He I per gram of materials: 



1 



- (riH + n e + TlHe)kTM + Ui H£iH + Up £H + Ui Hcl U Hcl + riHell CHcI 



(17) 



Here we have assumed that all He ill have already recombined to Hell, and so the total number density of particles is 
+ n e + riHc- £iH is the excitation energy of the i th excited state of H, eh is the ionization energy of H, tiHei is the excitation 
energy of the i th excited state of He I, and tHei is the ionization energy of He I. 
Recalling the definition of /h c = ^Hc/^h, we rewrite Eg. 1171 as 



§ (1+ X. + /He) — 

2 m H + /hoTOho 



Tm ■ 



+ 



/ho — Xi I 



Tli Hcl 



i 

3;HcII £HcI 



■ SiH 



Xp £H 



■ e. Hoi + ■ 



TO,H + /He^He t— 1 (/He — £iHcl)mi mH + /lie^He 

i 

The excitation energy terms in Eo. 1181 are found to be negligible 8 , and so we have 

6H 



de = |(l + a; e + /He) —7 

2 m-H + /EcmH 



+ 



k 



-T M + 



dT M + 



£HcI 



2 mH + /HeTOH, 
dxHcII ■ 



■T M + 



rnu + /Hcina 



dx v 



2 m H + fneinne m H + /H e m H e 

On the other hand, we can also write for an adiabatic process 



(18) 



(19) 



(20) 



7 By Friedmann equations, pc 2 a + 3(pc 2 +p)a = 0. Since the baryon pressure is much smaller than the energy density pc 2 , even though 
the pressure of baryon will be halved during the recombination process, it is still a good approximation to ignore the pressure. 

8 To be explained in the appendix. 
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where p is the pressure, given by the ideal gas law as 

k 

P = (1 + /He + X p + XHell) P -, Tm ■ (21) 

mn + /HeTlHe 

Now we may eliminate p to obtai n an equation with p, Tm, x p , XHcii onl y. Note that our method is similar to that in 
section 14 of lMihalas fc Mihalasl (11999!) . but there is an important difference. In lMihalas fc Mihalaa (I1998rl . the Saha equation 
is used to eliminate the ionization fraction. We certainly should not use the Saha equation in this case, but instead we will 
keep the derivative of the ionization fraction. Putting Ea. 12 II into Ea. 1201 to eliminate p, we get 



de 



(1 + /ho + x p + xnaii) — ; Tm 

mn + jH C mnc 



which can then be combined with Eg. 1191 and p cx (1 + z) 3 to obtain the modified temperature equation 
1 + z dT M _ 8a T U n e T M - T R / 2 e H \ l + z dx p 



^»-(-§&)r 



T M dz 3H(z)m e c n e + n H + nHe Tm V 3 kT M / 1 + /h c + x e dz 

(23) 



1 + - £HcI1 ^ 1+Z d^Hell 

3 fcT M / 1 + /ho + x e dz 

Eq. 1231 should then be used instead of the original equation Eq. Ilfil Hereafter we call the first term on the right to be 
the Compton Scattering term, the second the original adiabatic cooling term, the third and last the modified hydrogen and 
helium terms respectively. 



4 RESULTS 

Now we investigate the effects of recombination softening and reheating. The model that we consider is the ACDM model, 
with dark energy, CDM, baryon and radiation as the co mponents in the universe. Our choice of the cosmological parameters 
bases on the Running Spectral Index Model of WMAP jBennett et al 1 l2003l: ISpereel et ai]l2003h . In this model, the universe 
is flat, and so the total density parameter fio is 1.0 and Qk is 0.0. The contribution of the dark energy Qa is 0.73, and the 
total matter density S7m is 0.27, in which the CDM density Sic is 0.226 and the baryon density S1b is 0.044. The EOS index 
of dark energy is equal to —0.78 and is fixed. We use h = 0.71, To = 2.725 K and Y p = 0.24. The scalar spectral index n s is 
0.93, and the derivative of spectral index dn s /dlnA; is —0.031. We also simplify the calculation by assuming no reionization, 
lensing effect, tensor perturbations and massive neutrino in both models. Besides, we use COBE normalization throughout 
our calculation. 

We have modified the code RECFAST in CMBFAST 4.3 JSeliak fc Zaldarriagalll996h . The recombination process of 
He ill to Hell is calculated using the Saha equation. We solved a set of three differential equations, namely, the ionization 
equations for H and He I and the modified tempe rature equation for matter Eq. 1231 The equations are then integrated by 
using the integrator DVERK ifoeager et alJll99flfr . As in the original RECFAST, we turn off Hell whenever its ionization 
fraction is small enough. We use a fudge factor 1 in our calculation. However, we find that the conclusion that we make below 
is unchanged if we use the fudge facter used in RECFAST. 

Fig-0shows the values of the various terms on the r.h.s. of the modified temperature equation Eg. 1231 The recombination 
reheating tends to delay the recombination process. However, Compton scattering keeps the matter temperature close to the 
radiation temperature, because the ionization fraction is still quite large. This explains the shapes of the curves in Fig. Q 
This tight-coupling effect in fact has reduced the effect of recombination softening and reheating. 

The sum of all four terms in Eq. 1231 is plotted in the upper panel of Fig. [5] At large z, the coefficient is equal to 1. This 
can be understood by knowing that matter a nd photons co uple very well before recombination. Because the heat capacity of 
radiation is much larger than that of matter <|Peebles|ll993F) . the matter temperature will follow that of photons. The value of 
the coefficient is thus equal to 3(7 — 1) = 3(4/3 — 1) = 1. The coefficient deviates from 1 as z decreases. That is because the 
ionization fraction drops so fast that the Compton scattering can no longer force matter to follow the photon temperature. 
The percentage difference in (1 + z) d(hiTM)/dz with and without recombination softening and reheating effects is shown in 
the lower panel of Fig. [5] When the recombination starts, the sum of the terms is smaller than the original calculation; this 
means the temperature decreases more slowly after the modification. 

The percentage difference in the ionization fraction with and without recombination softening and reheating is plotted 
in Fig. El The original ionization history is shown as a subpanel. As the redshift decreases, the ionization fraction x e drops. 



9 Since the recombination time scale is long compared to that to establish thermal equilibrium in the cosmic plasma, the composition 
of the baryon gas remains nearly constant. Due to this fact and also the fact that kinetic energy of baryons during recombination is 
much smaller than the excitational energy of hydrogen, the ideal gas law is valid. However, the adiabatic index is changed because of 
the changing degrees of freedom during recombination. 
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Figure 1. The z dependence of the various terms of the modified temperature equation Eg. 1231 The solid, dashed, dot-dashed, and dotted 
lines are the Compton Scattering (CS), adiabatic cooling (AD), and the modified Hydrogen (H) and Helium (He) terms respectively. 




The terms dx p /dz and dxncii/dz in the new temperature equation are not zero for decreasing ionization fraction, and so the 
ionization fraction is also modified. 

With the ionization fraction at different redshift z, we can then calculate the visibility function g(z) = exp(— T)dr/drj(z), 
where r is the optical depth, r](z) is the conformal time. The peak of the visibility function defines the time of recombination, 
whereas its width defines the thickness of the last scattering surface. The visibility function and its changes are shown in 
Fig. [I] which clearly shows that the softening and reheating effect has delayed the recombination process. 

We can understand the effect of the delay in recombination in terms of entropy. The specific entropy s (entropy per unit 
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Figure 3. Change in the ionization fraction due to recombination softening and reheating. The subpanel shows the ionization history 
from between present and z = 10000. 
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Figure 4. The visibility function g and its change A3 when recombination softening and reheating are inluded. A fudge factor of 1 has 
been used in the calculation. 



mass) is defined by Tds = de +pd(l/p). For matter, the change of specific entropy in the recombination process is 



ds = — 



de- 4 dp 



fc(l + X p + IHcII + /He 



3 dT M _ dz 



3 £h 



2 T M 1 + 2 V2 fcT M / 1 + x p + xhcii + /a 



da; K 



(24) 



, 2 kT M ) l + x p + Xft e ii + /hc 
A certain amount of entropy is released when the ions and electrons recombined. The terms dx p and dsHcii correspond 
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Figure 5. CMB temperature and E-polarization spectra (upper panel), and the percentage differences of the spectra C; (lower panel), 
comparing results calculated with and without recombination softening and reheating effects. 



to this effect. The entropy is then shared among the particles remaining in thermal contact, as reflected in the prefactor 
1 + x p + XHcii + /h c . As the gas recombines, the entropy is shared among less particles, thus the gas is heated 10 and the 
recombination process is delayed. 

A delay of recombination changes the sound horizon at recombination, which directly affects the CMB power spectra. The 
new CMB power spectra and the percentage difference are shown in the lower panel of Fig. with the original CMB spectra 
shown in the upper panel for comparison. Since the recombination is only delayed slightly, the change of peak positions is 
negligible. However, the change of amplitude of Ci can be as high as 1% to 2%. In the results of WMAP, the temperature 
fluctuations At at the first peak , the next trough, and the second peak are 74.7 ± 0.5 fiK, 41.0 ± 0.5 /jK, and 48.8 ± 0.9 /jK 
respectively jBennett et alJl2003i) . They are equivalent to uncertainties of 1.3%, 2.5% and 3.5% in Ci at I = 220.1, 411.7 and 
546. Therefore, the error of calculation due to ignoring the recombination softening and reheating effect is about 10 times 
smaller than error of WMAP. However, WMAP is an ongoing project and the observational uncertainties will be reduced. 
Also, the next genera tion of CMB observations, such as PLANCK, may achieve precision better than 1% for I < 1000 
feersanelli et alJll99rl . and so the recombination softening and reheating effects should be included in future analysis. 

Besides experimental error in Ci, there is a statistical error know n as cosmic variance. T he theoretical limit of precision 
in the determination of Ci is ~1% at I = 100, and ~0.1% at / = 1000 JHu &: Dodelson|l2002t) . The effects of our modification 
are thus much larger than the cosmic variance for most I. Also, the effect of recombination softening and reheating increases 
as Qb increases. This is expected since our modification of the temperature equation is solely for the baryons. However, even 
when £Ib is as large as 1, the change in C; is only about twice as in the ACDM model. Thus our conclusion is the same in a 
large range of Ob- 

The change of the CDM transfer function T(k) is of order 0.01%. We normalize the transfer function at large scale to 
unity, so by definition there is no change in the large scale (small k). As the scale decreases, there is a region of T(k) with 
signature of fluctuations in photon-baryon fluid. The delay of recombination will change the sound horizon and so will modify 
the fluctuations. At even smaller scales (large k), the fluctuations are damped b y Silk damping. The damping length is about 
the geometric mean of the horizon size and the mean free path |Peacocj]|l99i) , and it is therefore larger in our calculation. 
This explains the decrease in the transfer function that we have found in our calculation. The new and old transfer functions 
approach a constant ratio because our calculation does not change the asymptotic behavior of the transfer function. The 
change in the transfer function is too small and can be safely neglected. 



One can also say that the cooling process is slowed down. 
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5 SUMMARY 

We have presented all results using a fudge factor of 1 in the recombination calculation. We have also repeated the calculation 
using the fudge factor used in RECFAST, and the conclusion is essentially the same. Whether the fudge factor itself should 
be changed when taking recombination softening and reheating into account awaits further calculation. 

The relative difference of the ionization fraction x e calculated by the modified RECFAST with the new temperature 
equation and the original RECFAST increases as z decreases for z < 300. However, the temperature equation breaks down 
for small redshift. As the universe continues to cool down, complications such as the formation of molecules, formation of 
structure at different scales, and reionization dominate the physics. 

In summary, we have studied the effects of recombination softening and reheating of the cosmic plasma on the ionization 
history, visibility function, the CMB spectra, and the CDM transfer function, using a generalized adiabatic index and standard 
cosmological parameters. With the standard cosmological parameters (ACDM model), the effects on the CMB spectra and 
ionization fraction are about 0.5% at I = 1500 and about 1% at I = 3000. These numbers are small compared to the 
current observational uncertainties but are not negligible. As the experimental errors become smaller, this improvement of 
the calculation should be used for determining the cosmological parameters. The effect on the CDM transfer function is much 
smaller, and so this modification is not important in the calculation of matter power spectrum. 

The work described in this paper was substantially supported by a grant from the Research Grants Council of the Hong 
Kong Special Administrative Region, China (Project No. 400803), and a Direct Grant (Project Code: 2060248) from the 
Chinese University of Hong Kong. 
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6 APPENDIX 

6.1 Comparsion of the energy terms 

For the density and temperature of the baryon gas at recombination, the major effect of line broadening is thermal broadening. 
The gas is nearly in thermal equilibrium and so follows Maxwell's distribution. The broadening is therefore about the width 
of the gaussian function. See lLansl Jl980ll for a more detailed discussion of the effect. 

The thermal broadening can be calculated by u(2kTM/muc 2 ) 1 ^ 2 , where v is the photon frequency. ISeaeer et alJ f2000) 
has shown that at recombination the energy separation between the n — 300 level and the continuum of hydrogen is nearly 
equal to the energy width of the thermal broadening. The electrons in level higher than 300 will then be ionized thermally. 
Simila rly we should inclu de only 300 levels in our calculation of the excitation energy of hydrogen. For helium, we should 
follow ISeaeer et all (l200rJl to use 200 levels. Without this truncation of states, the sum in the excitation terms will diverge 
unphysically. 

After the Hell recombination, the specific internal energy of baryonic matter is 
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Figure 6. Different terms in the specific internal energy. The lines represent contribution from thermal energy of all particles (TE), 
ionization energy of hydrogen (IE H), ionization energy of helium I (IE He), excitation energy of hydrogen (EE H), and the upper limit 
of the excitation energy of helium I (EE He). 



«H (m H + /HcTOHo 



71h(1 + Xp + XHell + fHe)kTM + ^] ^iHEiH + n p eu + niHcieiHcI + TlHcIieHoI 



(25) 



where the first term is the thermal energy of all particles, the second and the third terms are the excitation energy and 
ionization energy of hydrogen atom respectively, and the fourth and fifth terms are the excitation energy and ionization 
energy of He I. 

We are able to calculate the two ionization energy terms. We now simplify the hydrogen excitation energy term. When H 
starts to re combine, He has nearly re combined completely. The number density of H atoms with the electron in the i th state 
is given by jMihalas fc Mihalaslll999l) 

n m = 2Ci 2 n e n p T~ s/2 e ea/ * 2kTM , (26) 

with C = {l/2)(h 2 /2Tvm e k) [i/2 . 

The excitation energy in the i th state of hydrogen is given by e,H = (1 — l/i 2 )eH, and so the excitation energy term of 
hydrogen is 

1 \ ^ 1 — x p \ / h,h 



i . r \ / n »H£iH = ; — ;: > | EiH 

riH (rriH + /HcfriHc ) ^— ' m H + jEcmne '-^ \ «-h 

i i v 

= 2C Pcr » B (l - YJen i + 3 3/2 y 2 _ £H /^ fc T M (2?) 

m H + /HemHe TOH ^— ' 

i 

where p cr is the critical density. 

The excitation energy for He I is much more difficult to calculate. The helium atom is not hydrogenic. There is no close 
form for the energy states. One has to use methods like perturbation method or variational method. To calculate the energy 
states with high accuracy, one can refer to iDrakel l|l996l) . 

However for our purpose of checking the validity of neglecting the excitation terms, we do not have to calculate the 
correct energy states. Instead we find an upper limit of the states. We calculate the He I excitation energy with similar steps 
as the hydrogen. Then we construct a hydrogenic atom with a helium nucleus and one electron. Because we have neglected 
the repulsive force of another electron, the calculated energy in a certain state will be larger than the energy in the same state 
of the real helium atom. Therefore we can find a upper limit for the He I excitation energy. 

Fig. shows the five terms of the specific internal energy in Ea. 1251 The two excitation energy terms are much smaller 
than the thermal energy of all particles and the ionization energy of hydrogen, and so it is appropriate to neglect the excitation 
energy terms. The He I ionization energy term is kept because it is important before and during the He I recombination. 
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